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1.0 INTRODUCTION 

Knowledge of sea surface temperatures (SST's) is extremely useful in many 
applications including com mercual fishing, weather prediction and military uses. In 
the last two decades SST's have been determined fVcm satellite-based radiometers 
by a number of authors including Fritz and Winston (1962), Smith et al (1970)i 
Prabhakara and Dalu (197^)i McClaiTi and Abel (1977) and McClain (1930). Until 
recently, thermal inlVared measurements in the 10.5'* 12.5 /urn window region were 
mainly used for deriving SST's. In this spectral region the ocean surface behaves 
like an almost peii'ect blackbody radiator. The emissivity of the ocean has been 
estimated to be in the range of .967 to .997. By assuming it to be unity one can 
invert the Planck equation and compute SST. 

Two basic problems limit the accuracy of SST's derived ftx>m thermal IR 
measurements. One is the presence of clouds in the field of view (FOV). Clouds, 
being at higher altitudes than the ocean surface, absorb the radiation from the 
ocean surface and reemit radiation cotresponding to lower temperatures. Because 
of this even a small fractional covemge of the FOV will lead to unacceptable 
errors in the SST. It is, therefore, necessary to eliminate doud-contaminated 
FOV'sfom SST computations. 

The other problem limiting the accuracy of SST computations is the effect of 
the atmospheric constituents on the radiation emitted by the ocean. The 
atmosphere attenuates the radiation so that Ts, the sea surface temperature, is 
related to Ta, the equivalent temperature at the top of the atmosphere, by 

Ts = Ta dT (1) 

where AT is a correction term due to the atmosphere. This correction term is 
dependent mainly on atmospheric water vapor and could be greater than 8K in the 
tropical and subtropical are”^ Currently a regression procedure is used in the 
NOAA operational SST i»x>duction scheme GOSSTCOMP which uses the 10.5- 11.5 
^m measurements fx>m the Advanced Very High Resolution Infared Radiometer 
(AVHRR) on the TIROS-N series of NOAA satellites and the water vapor 
information fDm the High Resolution Ihfared Sounder (HIRS) for the atmospheric 
correction (Walton, et aL, 1976). This technique of finding the atmospheric 
correction has the drawback that the HIRS data provides relativdy crude estimates 


of temperature and humidity profiles and has a coarse spatial resolution compared 
with the AVHRR spatial resolution (McClain, I 98 O). 

It is, therefore, necessary to look into the possibility of obtaining the 
atmospheric correction bo the SST*s derived fVom the thermal IR measurements. 
MultispectrallR measurements have been proposed for obtaining this correction by 
Anding and Kauth (1970), McMUlin (1975), and Prabhakara et al (1974). In the 
present study a new method of obtaining the atmospheric correction is examined. 

Microwave measurements fhsm the Scanning Multichannel Microwave 
Radiometer (SMMR) on board the Nimbu9-7 satellite at^ used to derive the 
precipitable water which is then used to obtain the atmospheric correction for use 
with AVHRR thermal IR measurements to obtain SST's. The microwave 
measurements are minimally affected by the presence of clouds so that they 
provide all-weather capability for the determination of SST and precipitable water. 
Thus the atmospheric correction, although at a relatively coarse resolution, is 
available over cloudy as w^ as cloud flee areas. In this study, SST's derived fVom 
SMMR measurements are compared with SST's derived ftom both AVHRR and 
SMMR measurements, as well as with the operational AVHRR SST data. 


2.0 COMPUTATION OF ATMOSPHERIC CORRECTION 

2.1 Atmospheric EfYects 

As mentioned before, the temperatures observed by a sateUlte-based 
sensor are generally lower than the true surface temperature due to absorption and 
ree mission by the atmosphere of the surface radiation. The atmaopheric 
constituents causing this attenuation are mainly carbon dioxide, ozone, water vapor 
and aerosols. Experimental measurements by Burch (1970) show that in the window 
region absorption by CO 2 is negligible. The same is true for ozone also. For the 
10.5 ^m to 12.5 region absorption and ree mission by water vapor greatly exceed 
that by any other species except under very dry and hazy conditions (Blgnell et al, 
1963 ). The total water vapor absorption coefficient in the 10.5 wm - 12.5 Aim 
region can be written as (Prabhakara et al, 1974) 

K = Ke(v,T)+ Kp(v,T)+ Ki(v,T) (2) 

where and Kp are associated with the water vapor continuum and Ki with the 
water vapor lines. Experimental measurements by Burch (1970) and Bignell (1970) 
show that water vapor absorption depends strongly on the partial pressure of water 
vapor. Burch and BigneH have given the temperature variation of the mass 
absorption coefficients for water vapor. Their data can be used to compute the 
absorption by numerical integration of the equation of radiative transfer. This wiU 
enable us to compute the radiance transmitted £h>m the bottom to the top of the 
atmosphere. We can then define the atmospheric correction AT as the difference 
between the surface temperature Ts and the equivalent temperature T^ at the top 
of the atmosphere for a given slant path through a doud-ftree atmosphere. 

2.2 Equation of Transfer 

The basic equation of radiative transfer for an atmosphere in local 
thermodynamic equilibrium is 

-^ii-=Iv-Bv(T) (3) 

Kvpds 

wheire I v is the radiant energy in the ft^quency interval (v ,v+dv), Kv is the 
absorption coefficient, p is the density, s the distance variable and B the Planck 
function. The formal solution of equation (3) is (Chandrasekhar, I960) 


{ 


= Ioe-T(3»o)4.J' B(s,T)e-T(s^’) Rpds 


(H) 


where x (s^’) is the optical thickness of the material between the points s and s’ 
given by 


t(s^) = 




s’ 


Kpds 


(5) 


By using equation (5) and defining dm s pds the general solution to the equation of 
radiative transfer becomes 
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( 6 ) 


m is the water vapor mass along the path. 

Following Cogan et al (1974) the atmosphere may be assumed to be 
made up of suffLciently thin layers in which variations are linear. Then 


^ “ -^^-l ®*P (”^n Mp) + B(1 - exp(-kn Mp)) (7) 

where is the water vapor mass in a ccQumn of unit cross section in the layer n. 
k is the mass absorption coefficient in the 11.5 um window (Bignell, 1970) 
consisting of two parts k^i and k 2 < k^ is modified by changes in atmospheric 
pressure P ftx>m a standard value of 1000 mb, and k 2 is altered by changes in 
partial pressure e of water vapor fVom the standard value, k fcx* the n^ layer may 
be written 


kn=kil!L + k2fJL^ (8) 

1000 1000 

BigneQl (1970) found that absc:»?ption by foreign broadening of water vapor lines 
(with coefficient k^) has a positive temperature dependence of about 0.5% per 
and the e-type absorption (with coefficient k 2 ) has a strong negative dependence of 
about 2 % per °C. Hence equation (8) can be written as 

kn = f l(T) ki / Pn. V fad) kg/ flL- ^ 

V 1000 j V. 1000/ 


( 9 ) 


where f^d) 5 1 - .005 (303-Tn) (10) 

and f 2 (T) = 1 + .08 (30>Tn) (11) 

Tn is the mean bepperatio'e in the nth layer. 

2.3 Computation of W at er Vapor Correotion For The Thermal Window 

Channel 

It is possible to compute the radiance observed by a satellite-based 
sensor at the top of the atmosphere iyx}m equation (7) provided temperature and 
water vapor pressure are available as a function of atmospheric pressure fyx>m the 
ocean to the top of the atmosphere. From this radiance the equivalent blackbody 
temperature at the top of the atmosphere, (equation 1) can be obtained by 
Inverting the Planck equation. The correction T is then given by the difference 
between the surface temperature Ts and the temperature at the top of the 
atmosphere T^. 

In this study, the total atmospheric water vapor obtained fVx)m SMMR 
measurements are used for predicting the atmospheric ccxrection AT. 
Temperature and water vapor profiles are not available to evaluate this correction. 
Hence, it is necessary to determine the flmctional r^tionship between the 
correction aT and the total atmospheric water vapor W. Such a relationship can be 
obtained if a set of representative profiles of temperature and water vapor are 
available to compute both aT and W, For this study, due to the limited resources 
available, it was decided to use radiosonde data from ship soundings. These 
soundings were available on five tapes and covered the period February 11, 1979 
through March 17» 1979* Altogether there were 91 useable soundings covering a 
wide range of latitudes and longitudes. These profiles were used to compute both T 
and W for a number of slant angles, 0, by using equations (7 and 1). The 91 sets of 
at and W, as well as the temperature for the lowest point of the profile T^ were 
then input to a multiple regression routine which was used to determine the 
coefficients for a number of fmctional models. Eventually the following model 
was chosen on the basis of the multiple correlation coefficient and the sum of the 
squares of the residuals. 

A T = Sq + a*) • W + a 2 • Vj2 + S 3 W sec 0 
+ a^ • W • sec2 0 a5 T© 


( 12 ) 



The multiple correlation coefficient for this equation was found to be . 98 . The 
various coefficients in this equation and their standard deviations are given in 
Table I T@ in equation (12) is considered to be quite close to the surface 
temperature Ts. 

It is now possible to compute the atmospheric correction AT, given the 
precipitable water W and surface temperature Ts* 


TABLE I 


COEFFICIENTS OF EQUATION (12) AND THEIR STANDARD DEVIATIO NS 


Coefncients: 

a© 

ai 


®3 

34 

as 

Value: 

-4.315 

0.6666 

0.05648 

-0.2718 

-0.01603 

0.01582 

a : 

0.440 

0.0940 

0.00344 

0.1216 

0.04036 

0.00130 



3.0 COMPUTATION OF SEA SURFACE TEMPE^ T URSS FROM AVHRR AND 

SMMR MEASUREMENTS 

3.1 AVHRR Da ta 

The Advanced Very High Resolution Radiometer (AVHRR) is a cross- 
track scanning radiometer with the spectral channels shown in Table E. The IFOV 
of the sensor is approximately 1.0 km at the sub-satellite point. For this study the 
Global Area Coverage (G A C) data at the reduced resdutioh of M km was used. This 
data was used in two forma. The first was the raw GAG 4 km data» from which, for 
each pixel, the latitude, longitude, solar zenith angle, an^ of observation through 
the atmosphere, bdghtneases fhom the visible channel, and the unoorrected 11 
micron IR temperature can be figured. The temperature corrections are done 
starting with this IR data. 

The other form is that of the operational NO A A GOSSTCOMP SST data 
derived fhom the AVHRR and HIRS data. This SST data, which is retrieved fhom 
observations averaged over 50 km target areas and stored on tapes covering seven 
day periods, is used in this study only for comparison with the corrected 
temperatures obtained fhom the raw AVHRR data. 

3.2 SMMR Data 

The Scanning MultichanndL Microwave Radiometer (SMMR) (Figure 1) 
provides orthogonally polarized antenna temperature measurements at five 
microwave fl?equencies, 6.6 GHz, 10.7 GHz, 18.0 GHz, 21.0 GHz and 37.0 GHz with 
an absolute accuracy of leas than 29K (rms) at each fV^quency. The antenna's main 
beam is offset 42° fV^m nadir and scans the earth in a conical pattern with a half 
ang^ of about 25^ (Figure 2). The SM M R data used for this ^^dy was fhom SMMR 
CELL-ALL tapes, a product of the second stage of SM M R processing. Each logical 
record contains earth-located SMMR brightness temperatures for a block 780 km 
by 780 km on the surface of the earth. The antenna temperature measurements 
fVx>m the SMMR are converted to brightness temperatures and averaged into grids 
of various sizes based on the fV^quency of the observation. Table HI shows the grid 
size, the number of bands and cells in each grid, and the frequencies of the 
observations which api^ for each grid. 

The SMMR brightness temperature for each c^ of the 780 km x 780 km 
block have been corrected for the effects of the antenna sidelobe pattern. In 
addition the CELL-ALL tape brightness temperatures have been corrected for the 



TABLE H 


THE AVHRR CHAMNEL SPECTRAL RANGES 


CHANNEL 

SPECTRAL RANGE 

REGION 

1 

0.55 - 0.68 

. VisLUle 

2 

0.725 - 1.10 

VMhLe/ N ear IhfV'ared 

3 

3.55 - 3.93 

Ihll^ared 


10.50 - 11.50 

InlV'ared 


SMMR INSTRUMENT CONFIGURATION 


a#- 



NIMBUS-G 

SMMR SO^N GEOMETRY 



I SWATH ANGLE, deg. -50 
SWATH Width, km - 778 



TABLE HI 


i 

FREQUE N C Y AND CELL SIZE OF OBSER V A TIO NS ON THESMMR CELL^ALL TAPE 


GRID 

GRID SIZE (km) 

BANDS X CELLS 

CHANNELS INCLUDED IN GRID (GHz) 

1 

156.0 

5x5 

37.0, 21.0, 18.0, 10.7, 6.6 

2 

97.5 

8x8 

37.0, 21.0, 18.0, 10.7 

3 

60.0 

13x 13 

37.0, 21.0, 18.0 

4 

30.0 

26 x26 

37.0 


0 


effects of polarization mixing between the vertically and horizontally polarized 
components of the signal received by the SM M R antenna. 

To retrieve atmospheric water vapor and (SMMR) sea surface 
temperatures, the algorithms developed by Wilheit and Chang (I 9 BO) may be used 
on the cell data. The water vapor algorithm uses the 37.0, 21.0 and 18.0 GHz 
channel information, which 1s avaUablfi for the 13x13 grid whose ceU size is 60.0 
km on each side as wdi as for coarser grids, while <the SST algorithm uses 
information fVx>m all channels, which is only available ihom the 5x5 156.0 km 
square gid. The 13 x13 60 kmSMMR grid was chosen as the reference grid, into 
which A VH R R pix^ would be registered, and SM M R SST's interpcQated. 

Only data fh>m the ascending (daytime) nodes of SMMR were used, so 
as to allow the closest time coincidence with the A VH R R ascending (daytime) data. 
The local time for an ascending SMMR node is always about noon (12:00 houre), and 
the local time for an ascending AVHRR node is always about 3 P*m. (15:00 hours). 
It was decided for this study that the resulting 3 hours difference is insignificant 
for 60 km square areas for those parameta^ measured by SM M R . 

3.3 Overview of Computation of SST^s 

Figure 3 gives a gross schematic of the computation of SST’s by using 
both AVHRR and SM M R data. First, total atmospheric water vapor W is computed 
n?om SMMR brightness temperatures. The precipltable wata? is then used to 
compute atmospheric correction as expilained in Section 2.0. The correction uses 
the brightness temperature T/^ for the 11.5 Am channel over dear areas lW>m the 
AVHRR GAC data (Section 3*5 regarding doud filtaing). T^ has not been 
corrected for atmospheric effects. It is used as a first approximation for the 
surface temperature. The aT so derived is then used to ccxrect T^ and obtain the 
first approximation of the carected surface temperature Ts. This is again input 
Into equation (12) and aT is computed a second time. This aT is then used to 
correct Ta again to obtain a final surface temperature (SST a) fVx)m AVHRR 
brightness tem peratures. 

Sea surface temperatures are also computed tirom SMMR txightness 
temperatures directly by using the algorithm of Wilheit and Chang (i 960 ) (Section 
3 . 2 ). These are denoted by SSTs. Finally, the coresponding sea surface 
temperatures iVom the NO A A operational product SSTq are extracted fVx>m the 
seven day GOSSTCO MF tapes. The three sets of SST's are then compared and the 
statistics of the residuals are computed. 



FIGURE 3 

OVERALL SCHEMATIC OF COMPUTATIONS AND COMPARISONS 
OF SST'8 FROM AVHRR AMD SMMR DATA 




3.4 Area Selection and Data Handling 

Six open ocean areas were initially selected, on the basis of the 20° x 
2Q0 open ocean blocks used in Kogut and Hopkins (I 98 O). SMMR CELL-ALL data 
covering the time flnora Febnaary 15 to March 17» 1979 were used to plot the 
approximate locations of the centers of aU 780 x 780 km SM M R blocks lYom each 
of the days covered by the tape, 

A VH R R visible and infhared imagery was then used to ascertain, out of 
all the locations of SMMR blocks within the initially selected areas, which 
locations were relatively dear. A little patchy doudiness was acceptable, so long 
as there were sufYldent areas to ensure that at least isolated pixels of infY^ared 
data would reach from the surface over most of where the 780 x 780 km SMMR 
block covered. Viside and IR GOES prints were used to supplement this 
information where possible. 

Baaed on this information, AVHRR GAG data was ordered for ten 
coverages of a region on the earth on a given day. After the sdected-coverage 
GAG data arrived, both AVHRR visible histogram runs (see Section 3.5) and 
temperature estimation runs (following the algorithm of Section 3.3) were made. 
In each run, the SM MR tape was scanned first, until an ascending>node (daytime) 
block within the given region and day was found. An array re{H:'esenting a grid of 
latitude and longitude lines covering the 780 x 780 km SM MR block was then filled 
in with information about how far over and how far up on the SMMR block each 
grid point was. This was done byr 

1. Figuring what the proper* rotation and stretching algorithms 
would be, were the earth's surface flat, fVum the 
latitude/longitude grid to the mesh of 60 x 60 km cell 
centers in the SMMR block, based on the latitudes and/or 
longitudes of three cell centers. 

2. For each cell, obtaining the residual between its true 
latitude/longitude and that computed fnom the inverse of 
the transformation fhom (1) above. 

3 . For each point in the latitude/longitude grid, 

(a) finding its supposed SMMR block portion Ororn 
the transformation of (1); 



(b) using that position to interpolate what residual 
there would be ftx)ra the inverse transformation 
fhom around that position to latitude and 
longitude; and 

(c) adding that residual (calculated minus real 
latitude/Qongitude) to the latitude/longitude of 
the latitude/Longitude grid point and obtaining (a 
reasonably) true SMMR block location by using 
the (forward) transformation of (1). 

In temperature runs, the SMMR water vapor values were computed 
directly fVom the 60 x 60 km cdl data, while the sea surface temperatures were 
computed fh)m the 156 x 156 km data. A polynomial was then fitted over the 
SST's, and was then evaluated at the center of each 60 x 60 km 

The AVHRR tape file for the corresponding region was then read, and 
for each AVHRR pixel, the corresponding 60 x 60 km SMMR ceU was found by 
interpolating the pixd.'s latitude and longitude on the grid for SMMR block 
location. In histogram runs, the visible values were used if the pixel was within the 
780 X 760 km block at alL In temperature runs, the IR value was added in as part 
of an averaging process for every 60 x 60 km celL For all runs, the final 
calculations and printouts were done after the AVHRR tape file had been 
completeOy read. 

3.5 Cloud Filtering of AVH R R Data 

The identification of cloudy pixels was done oy establishing brightness 
thresholds for two of the ^dsibLe "reflectance channels^'. The two channels used 
were the red (0.725-1.10 pm) and green (0.55-0.68 pm) of the AVHRR. The 
radiances in these channels were normalized by dividing by the cosine of the solar 
zenith an^. 

Special runs were made to obtain histograms of the radiances in each 
channd, which were then plotted on the printer output. A typical pair of 
histograms are shown in Figures 4 and 5. 

The large peak at the low radiance end of the histogram cotresponds to 
reflection fhom the oceaui surface. Peaks at higher radiances correspond to 
reflection fTom cloud-ocean combinations or from clouds. The cloud-no cloud 
thre:^QLd was sdected to be Ro -t- AR where Ro is the radiance corresponding to 


iUSTOGRAM OP YB>IIED 




• • • HI • 0 9 < 


Sill ISSIS 


•^S 7 Si tr« SCR K}« RilQ 73 


HOOO OOOOOOeOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOQOH 


8S8 8e99S9oaoeeeCtoooooooooooooooeooeooooooooooooe 

- » 5vrf ^;»^ 53 ;; 5 - }>J 5» 55 » » SS 




0 

0 


TXnCAL USTOCRAN OP TRB RADURCESn THE RED (0.72>1.10 UB) RVHRR CHANNEL 


IHBTOGRAM OP nS-GREEN 







$ ; ^ ^ K2HQ s ai 73 ^9 79 


HOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOH 

1 


888889999 Oooooeoooeeoeeoeeooooooeoeooeoooooooooo 

^ fi Sg SJ tS ffSS a JCa«« JPg&S 


000 


eo 


£::::: SSSii!ISIsiSSsS 33 S§S£:SS 3 :§§»Ss»»ss!:i:$!s 3 aa 5 




nCURE 5 

TmCRL RBTOCRiM OF THE RASUMCESIR THE CHEER (0.5^4^ ub) ATHRR CHARREL 


the peak representing reflections ftx>m the ocean and A R is determined fVx>m the 
ft^quency distribution between this peak and the cutoff of the lower radiance half 
of the histogram. AR was chosen conservatively so as to preclude even small 
fractions of doud coverage in the FOVs. 

The threshold for each visible channel was then entered as input into 
the temperature runs, which would compare the normalized radiance fhom each 
channel for each AVHRR pixel with the corresponding threshold, and use the IR 
fVx>m the pixel only if each radiance was less than the corresponding threshold. 


4.0 RESULTS 

Seven SM M R blocks of 780 x 780 km have been used for analysis in this study 
(Figure 6). Tables IV(a) through IVCj) show some of the results for the block 
centered at 279.9°S laUtude, I 8 . 9 OW longitude, on February 27, 1979. The 13 x 13 
eieraents of Tables IV(b) through IVCj) refer to the SMM R 60 x 60 km ceQls over 
which the average water vapor is retrieved. Table IV(a) shows the SST’s retrieved 
from SMMR brightness temperature data over the 5 x 5 km grid. Table IV(b) 
shows the interpolation of this over the 13 x 13 60km grid. The temperatures seem 
to be very low for latitudes in the order of 28^S. Table IV(c) shows the total water 
vapor for each 60 km cell. The values range ftxjm 2.7 gm/cm2 to 4.7 gm/cm^. 
These water vapor values are available even for those pixels which are defined as 
cloudy by visible threshold. Table IV (d) shows how many AVH RR pixels were used 
for each 60 km celL These numbers show where there are areas of partial 
cloudiness and complete overcast. 

Table IV (e) shows the brightness temperature of the A VH R R thermal channd. 
fcx' each i^eL These values are, as yet, not corrected for atmospheric water 
vapor. They are on the order of 29 1®K. Table IV (f) shows the final corrected SST 
computed by adding the atmospheric correction derived fmm SMMR data to the 
brightness temperatures derived iVxcm AVHRR data. It may be observed from 
Table IV(f) that the SST values are on the order of 296^K, considerably higher than 
the SMMR derived SST*s. Table IV(g) shows the SST's lh>m the NO A A operational 
product tape which were gridded on the SM M R 60 km grid by the nearest neighbor 
approach. It may be seen fVom Table IV(g) that the values are dose to the SST 
derived ihom AVHRR and SMMR. Tables IV(h) through IV(J} summarize residuals 
between different methods of deriving the SST, as follows: 

TablelV(h): SST(AVHRR & SMMR)-SMMR(AVHRR OPERAnONAL) 

Tableiva): SST(AVHRR & SMMR)-SST(SMMR) 

TablelVO): SST(SMMR)-SST(AVHRR OPERATIONAL) 

It may be observed that the SST derived fhom AVHRR and SMMR is at most 
about a degree less, on the average, than the N 0 A A operational AVHRR, while the 
SMMR SST is about 21 or 22 degrees less than both of these. This, in combination 
with the latitude, would seem to indicate that the A VH R R-^M M R derived ^duct 
corresponds at least reasonably weU to surface truth, while the SMMR product is at 



LOCATimS OF THE 780x780 SMMR BLOCKS USED FOR ANALYSIS 




TABLE IV(a) SMMR SST'a (®K) RETRIEVED OVER. THE 5x5 GRID 


273.12 

276.57 

274.29 

277.37 

278.47 

274.26 

275.64 

275.10 

277.56 

278.98 

274.26 

273.34 

276.95 

275.35 

277.64 

274.56 

274.15 

275.14 

275.13 

277.45 

274.07 

273.55 

273.69 

274.31 

275.23 


TABLE IV(b) SMMR SS’nsCOE) INTERPOLATED OVER THE 13x13 GRID 


cA<T>C7>^oocoeo^'t^t*-'vp ir>u^ 

Nw» Nw 


f>^cr^(n(^3(^iCV(^j^3ro^^p^M(^J 


oo^-^^^^v^o^-o^c^l/^c^^oo 

NfviKjhjKjNKjN<\jKjNNKJ 


CT>oinc^rnr";sjt>«^Qw:\p«iap 

S5^5ll2^|§8^|is| 

rnrnrnpocvtvrnrOrornfnfnrsi 


r; q OJ «j g oq q in q ar. 


in <\j o g> u> q <\i J2J Vg jn 2 

o\o ovo®oo®2iq® 

\oCofnooq*^*"oSr-f?3cyvoo^ 

cr^MMnrS w c3Mc3Mr*So3o^nS 


^ oj(N^ qq nrjq q q 


^^^r■(n 


jTrntnrnrofnrniCsjrnpn 


2 rn (T> 
^ o q »- 

(nf»S(nr>S 


^.^.vo £T (M o>vo <M oo mao 
• •••••••••••• 

v^\ovo'XJ‘o*^voininin«riTm 


jT(nmMc3mmmmo'Smn^m 


r^qo^qr^qqqo^qqrjq 

vo^%ovo^oininininirjs-J3’m 

KiRjRiRiRiKiKifiKiRifiKiRi 


JT jqr m ol cvHM M r»^ ff^ M m M 


9 ft ^ 210 

rocp 


N Or-<UQ’“ 
m »“ r" r- m Q *“ 
oj qq m 
^ (»S M cj m M 


o^ ina- mm 
mmmmmm 


r- (\j m m m nj »- o s- in »- oo a- 
inintnininuSinin^jr^mcy 

Kl Jm ivT ^ ^ ?v ^ 5 m ^ 


inqq»^q»^qqooj^qq 

pr^mmrfSmojmrncr^o'SMmm 


jrt*-ooa>oooqt-vomot^m 
c\j::r^^ rncvj 


m 

q<M,r-U>t^qJT^ 9> (M 3 VO 

min^ fnSq^ n? KJ <\/ ^ q q fS 
oS (tS m r»*Kvj (V tvSf»^n^mmmm 


o mvovo t^ooco VO j» <v o CM 

> 

.< 

Sm^ 

oojJ9'in»^vDvdm^<Mo 

• • • • _f‘' • • • » _ g • 

JTflrj^-^^/^y^r-sirsr-^rarrnrvj 

RiSiRKiRlRlflSlRRlKiRlS 


3 T- 3:00 min^-ooocMf-j^o 
oom®pin.^'=rs'=9‘in»~fn 
qqqinc5qmmmi»’.»mr- 

j^.5rmMm(MMmo^(YSTvSn*io^ 


^t^qojmqqqqqo^qqq 
r- mar Jirarararararo-aro- m<M 

;:5MfM5M5M5M?M5M5M5M5M5M^ 


Yy^r-0>arqOVOOOO>*-r-fMvoirt«-. 
r-0^*J^■^ar^•N•'00r-^S^-»- 
^qqq^qqqmqinqmq 
t-aro’ar mmmmmmmmmiri 


qqqcMqqqqqmqqt^ 
cnrn^^^sTsri^iT:^sT mm 
rM(MCMCMC4(MCM<MCM?M<M<\JCM 


ini-ocMfMovrvjirr-t^oopmt- 

5-»“inirmr-v£>arcr»ocrvin;3-cM 

^»-tn{\jaotnmmarvDt^c^;3-CM 

C/J ••••••••••••• 

ST ST mmmmmmmmmm 


TABLE W(d) NUMBER OF AVHRR PIXELS USED IN EACH 60 KM CELL 
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TABLE WO) FINAL CORRECTED SST^FROK SMMR AND AVHRR MEASUREMENTS 
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TABLE IV(h) RESIDUALS BETWEEN SST(AVHRR 4 SMMR) AND SST(AVHRR OPERATIONAL) 
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TABLE IVCi) RESIDUALS BETWEEN SSTC5MMR) AND SST(AVHRR OPERATIONAL) 
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least 15 degrees low. This discrepancy in the SM MR data corre-sponds with the 
amount of bias in the SM M R cell tape data which others have noted. 

Table V shows the average results for all seven 780 x 780 km SM M R blocks. 
The average of the average residuals is -.149» and the average standard deviation is 
1.022. It may be seen that most values derived fVx>m the AVH R R and SM M R data 
coincide reasonably wdQ. with the N 0 A A operational product tape data. 
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5.0 CONCLUSIONS 

In this study, microwave measurements fH>m the Nimbus-7 SMMR were used 
to derive the atmospheric precipitabile water, which was then used to obtain the 
atmospheric correction for use with AVHRR thermal IR measurements to obtain 
sea surface temperature. The resulting SST’s compared were with the NO A A 
operational sea surface temperature measurements, and the two sets of 
measurements were found to be in reasonable agreement. For the limited sample 
analyzed in this study the average residuals between the two sets of measurements 
was .15^K with the N 0 A A operational SST’s being :^htiy greater. 

This demonstrates that the use of microwave measurements is a promising 
tool for correcting IR sea surface temperatures for water vapor effects even under 
cloudy conditions. This can be done either by using Nimbu»>7 SM M R and TIR OS-N 
or NOAA-6 or NOAA-7 AVHRR measurements and matching them, or by using 
instruments ftem any future satellites where both the scanning IR and microwave 
instruments are on board. To further check the accuracy of this algorithm, the 
results fVem it should be compared with actual surface truth such as that fVem 
buoys. Also, since the method of doud thresholding "by hand" used in this study is 
impractical for operational use, the effects on this sQgorithm of automating the 
doud thresholding procedure should be studied. 



